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1. INTRODUCTION 

Under the application of a sufficiently strong 
electric Field, a suspension may undergo reversible 
phase transitions from a homogeneous random 
arrangement of particles into a variety of ordered 
aggregation patterns. The surprising fact about electric- 
field driven phase transitions is that the aggregation 
patterns, that are observed in very diverse systems of 
colloids, display a number of common structural 
features and modes of evolution thereby implying that 
a universal mechanism may exist to account for these 
phenomena. It is now generally believed that this 
mechanism emanates from the presence of the long- 
range anisotropic interactions between colloidal 
particles due to their polarization in an applied field. 
But, in spite of numerous applications of the electric- 
field-driven phenomena in biotechnology, separation, 
electrorheology, materials engineering, micro-devices 
for chemical analysis, etc. which have expanded 
rapidly over the last decade, our understanding of these 
phenomena is far from complete. Thus, it is the 
purpose of the proposed research to develop a theory 
and then test experimentally, under normal- and low- 
gravity conditions, the accuracy of the theoretical 
predictions regarding the effect of the synergism of the 
intcrparticlc electric and hydrodynamic interactions on 
the phase diagram of a suspension. 

The main results from our theoretical studies 
performed to-date enable one to trace how the 
variations of the electrical properties of the constituent 
materials influence the topology of the suspension 
phase diagram and then, by using an appropriate phase 
diagram, to evaluate how the electric- field-induced 
transformations will depend on the frequency and the 
strength of the applied field and the particle 
concentration. 

2. ELECTRIC-FIELD-INDUCED PHASE 
TRANSITIONS 

The first step in the development of a theory for the 
electric-field-induced phase transitions in a suspension 
of electrically uncharged conducting particles 
dispersed in a conducting fluid was taken in Refs, fl, 
2]. To begin with, the density of the free energy of a 
suspension containing randomly arranged hard spheres 
subject to an electric field was written fl] as 
F = k B Tf 0 (c)/v p -W (1) 


where the first term in Eq. (1) refers to the free energy 
of a suspension in the absence of the electric field, as a 
function of the temperature T with k B being 
Boltzmann’s constant, and the volume concentration of 
the particles C with V p being the particle volume; the 

second term in Eq. (1) refers to the electric energy of a 
suspension caused by the interaction of the particles 
with the applied electric field and by the electric-field- 
induced long-range interparticle interactions. The 
concentration dependence of the function f 0 (c) in Eq. 

(1) is well understood [3, 4]. 

The main obstacle which had to be overcome in 
constructing a theory of electric-field-induced phase 
transitions was that, for materials such as conducting 
suspensions whose complex dielectric permittivity 

8* — e' — ie'' varies strongly with the frequency of 
the applied electric field, the electric energy density 
W in Eq. (1) cannot be constructed using macroscopic 
electrodynamics [5]. The reason for this difficulty is 
that the frequency dependence of the dielectric 
permittivity implies that such a material contains 
mobile charges and electric dipoles capable of 
orientation, so that the stored electric energy will 
depend on the time-history of how the electric field was 
established [6]. 

To be sure, an expression for the electric energy 
density W of a conducting material, namely, 
Brillouin’s formula (sec Eq. (2) below), can be derived 
from macroscopic electrodynamics [5, 6], but only 
when the energy dissipation in this material is 
negligibly small and the time-variations of the applied 
electric field are very slow compared to the rate of 
relaxation phenomena, i.e. 

(2) 

only if |e"((d)| « e'(cd) and ©I* <<: 1 where co 

and E((o) are the frequency and the Fourier amplitude 
of the applied electric field, respectively; the value of 
t s determines the relaxation time of dielectric 

phenomena; and ( ) denotes the time average. 

Recall, that the real component £’ gives the dielectric 
constant of a material whereas the imaginary 
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component z” determines the power dissipation (loss) 

in this material due to its conductivity. As can be seen 
from Eq. (2), Brillouin’s formula is applicable only 
under the severe limitations that effects of conductivity 
are negligibly small so as not to cause losses in the 
stored electric energy. That is why Brillouin’s formula 
cannot be utilized in Eq. (1) to describe electric -field- 
driven phase transformations in conducting colloids 
which strongly depend on conductivity effects. 

To overcome the limitations of Brillouin’s formula, 
we developed [1] a microscopic theory for the electric 
energy density W of conducting spheres dispersed in 
a conducting fluid when the particles were arranged 
randomly and provided that the particles and the 
suspending fluid can be described by the model of a 
leaky dielectric; i.e. when their dielectric constants, E p 

and £ f , and conductivities, G p and O f , are 

frequency independent. This model corresponds to the 
classical mechanism of the so-called Maxwell-Wagner 
inlerfacial polarization typical of colloids [4, 7]. In this 
case, the short-term polarization of the particles and of 
the fluid is determined solely by their instantaneous 
polarization whereas their long-term polarization arises 
from the build-up of charge at the interface between the 
particles and the surrounding fluid. We exploited two 
microstruclure-based techniques for solving the 
problem. One of them is based on a mean-field 
approximation (a cell model) where the average field 
acting on each particle is considered to be the well- 
known Lorentz-Lorenz local field [5, 6] rather than the 
applied field, while the second made use of statistical 
methods and a renormalization technique for 

calculating the 0(c 2 ) term the expansion for the 

electric energy density W in powers of the particle 
concentration. But, as we found [1], the main 
contribution to the electric energy of a suspension as 
well as to its complex permittivity is given by 
multiparticle interactions leading to a change in the 
local electric field acting on a particle. That is the 
reason why, for example, the cell model, which 
accounts only for this effect, yields the well-known 
Maxwell-Wagner expression for the complex 


permittivity of a suspension which correlates well with 
experimental data. 

The theory of Ref. [11 relates the electric energy 
density of a conducting suspension to the dielectric 
constants and conductivities of the particles and of the 
suspending fluid, the particle concentration, and the 
frequency and the strength of the applied electric field. 
To illustrate the advantages of this theory, the first two 
terms in the expansion of the electric energy density 
w in powers of the particle concentration were 
substituted in Eq. (1) which was then used to study the 
phase separation of a dilute conducting suspension 
subject to strong electric fields [11. On this basis, we 
then calculated in Refs. [1,2] the critical conditions 
beyond which the random arrangement of the particles 
becomes unstable in the presence of spatially uniform 
and non-uniform electric fields. 

The main objective of the present work is to extend 
our previous studies [1,2] beyond the dilute regime 
and to develop a microscopic theory for phase 
diagrams of concentrated conducting suspensions 
subject to strong dc and ac electric fields. To this end, 
we employ in Eq. (1) the full expression for the free 
energy density of a conducting suspension, as a 
function of a concentration, which we derived in Ref. 
[1] using the mean field approximation, and then 
investigate the topology of the suspension phase 
diagram, i.e. “the particle concentration-the electric 
field strength”, and relate it to the electrical properties 
of the constituent materials. 

As was shown in Ref. [\], the expression for the 
electric energy density W of a conducting suspension 
being subjected to a dc electric field E for short 
(t « t s ) and for long (t » t s ) times, is given by 
Eqs. (3) and (4), respectively, while, for an ac electric 
field E 0 COSCOt and for the long-term t » t s 
regime, the equation for the time average of the electric 
energy of a suspension is given by Eq. (5) 


W = 


e f I + 2cft E 
2 l-cP E 


(3) 


W = 


W = 
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For brevity, we included the vacuum permittivity £ 0 
as a multiplier in the definition of dielectric constants 
8 p and £ f in Eqs. (4)-(5). 

Equation (3) and the first terms in Eqs. (4) and (5), 
which are determined solely by the instantaneous 
polarization of the particles and the fluid, are consistent 
with the relationship for the electric energy of a non- 
conducting suspension W = “£ s (c)E 2 given by 

macroscopic electrodynamics [5] with £ s (c) given by 
Maxwell’s equation [7]. The second terms in Eqs. (4) 
and (5) correspond to the energy required to build the 
charge at the particle surface and to redistribute the 
electric field inside the particles and the suspending 
fluid when the time constants of the particles and of the 
fluid arc different; that is when £ f / CJ f ^ £ p /o p [1]. 


P =— — . 

’ Pa a p +2a f 

Consequently, as seen from Eq. (6), as long as the third 
term on the right-side is positive, in the other words as 
long as the interparticle interactions increase the 
suspension energy, an increase in the strength of the 
electric field decreases the value of 3ll/8c and 
finally renders the random arrangement of the particles 
unstable when 

M = o (7) 

dc 

Now, Eq. (7) represents the spinodal curve in an 
equilibrium phase diagram of particle concentration vs. 
the electric field strength of a conducting suspension 
subject to an electric field. Specifically, this curve is 
the locus of points for which the curvature of the 
suspension free energy changes from convex to 
concave [8], 

For the case when 3 2 W/3c 2 > 0 over the entire 


But since the second term in Eq. (5) approaches zero as 
COt s — > oo the electric energy of the suspension for 
high frequencies becomes the same as that for a non- 
conducting material. On the other hand, for a slowly 
varying electric field, 0)1 s « 1 , Eq. (5) yields the 
expression for the time average of the electric energy 
which is consistent with the expression given by the 
substitution of the Maxwell-Wagner relation for 

£'(co) into Brillouin’s formula, Eq. (2), [1]. 

Now, the osmotic pressure of a suspension Yl 
(equal to -(3F/0 V) nte [8]) and the chemical 

potential of a particle J 1 (equal to (F + FIV)/N [8]) 
can be evaluated from the equation of the free energy, 
Eq. (1). On calculating these expressions, f 0 ( c ) is 
conveniently expressed in the terms of the suspension 
compressibility factor Z(c) where a good 

approximation for Z(c) of a suspension in the 
disordered state is given, for example [3, 4], by the 
Camahan-Starling equation and the asymptotics 
Z— >°° as C — > C m where C m ~ 0.63 — 0.64 

corresponds to random close packing of spheres. 

The random arrangement of the particles in a 
suspension is stable as long as the osmotic pressure 
increases with concentration, which can be wntten as 


a 

3c 






>0 ( 6 ) 


range 0 < C < C m , the spinodal points exist over the 

entire concentration range as well, so that increasing 
the strength of an electric field applied to this 
suspension will eventually cause the electric- field- 
induced transition. Thus, the phase diagram of such a 
suspension consists of the low- field one-phase region 
which includes the random spatial arrangement of the 
particles in the absence of an electric field and the 
high-field two-phase region corresponding to the 
appearance of aggregates caused by the action of an 
applied field. This phase diagram appears to be similar 
to the phase diagram concentration vs. temperature of a 
binary fluid or a binary alloy with a miscibility gap [81, 
so that these species exist in solution at high 
temperature but their mixture eventually separates into 
coexisting phases below the critical point of 

miscibility. In this regard, the application of an electric 
field to a conducting suspension is equivalent to a 
quench of an atomic system from its high -temperature 
one-phase state. Although such a similarity between the 
structural ordering in colloids and the structural 
behavior in conventional atomic systems- gas, liquid, 
crystal, and glasses- under suitable conditions has 
already been well recorded [3, 4], electric-field driven 
phase transitions in colloids have not been treated from 
this point of view thus far. 

The critical point on the spinodal curve of the 
suspension, namely, C CT and E 2 r , coincides with the 

inflection point of FI as a function of C , so that Eq. 
(7) has a multiple root. In a two-phase region (at 
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E 2 > E^ ), the value of the osmotic pressure and of 

the chemical potential of the particles in the coexisting 
phases are equal. The latter leads to the following 
relations [8] between the particle concentrations in the 
both phases, C { and c 2 , respectively: 

H(c,)=^l(c 2 ) and 11(0,)= Il(c 2 ) ( 8 ) 

Now, Eqs. (8) represents the coexistence curve in the 
phase diagram of a suspension. This curve is located to 
the left of the spinodal curve in the low-concentration 
part of the phase diagram and to the right of the 
spinodal curve in the high-concentration part of this 
diagram. Moreover, C { — > C cr and C 2 — > C cr as 


On the other hand, if the third term on the right-side 
of Eq. (6) is negative, in other words when the long- 
range interparticle interactions lower the suspension 

energy (i.e. when 3 2 W/8c 2 < 0 over the entire 


range 0 < C < C m ), an increase in the strength of the 
applied electric field increases the value of 5n/5c . In 


this case, the random arrangement of the particles 
appears to be stable, so that the phase diagram of such 
a suspension reduces to the one-phase domain. 

As can be seen from Eqs (4)-(5), the sign of 

5 2 W/3c 2 depends on the particlc-to- fluid ratios of 


the conductivities and of the dielectric constants, the 
frequency of the applied electric field, and the particle 
concentration. For a high-frequency electric field when 
COt s — > co,Eq. (5) yields 

3 2 W 3e f p 2 E l 


dc 1 


2(l -P e c) 3 


>0 


(9) 


Thus, the sign of 3 2 W/3c 2 at P E 0 always 


becomes positive over 0<C<C m when the 

frequency of the applied electric field becomes 
sufficiently high. Hence it follows that there always 
exists a threshold value of the frequency, above which 
the particles having a dielectric constant different from 
that of the suspending fluid (regardless of the mismatch 
of their conductivities, O py /a f ) will aggregate as the 

strength of an applied electric field becomes 
sufficiently large. However, as we shall see, the sign of 

a 2 w/a c 2 may become negative for low frequencies. 
If this happens over the entire concentration range, 
0 < C < C m , the electric-field-induced aggregation of 

the particles will not occur at all. A peculiarity of 
concentrated suspensions is that the sign of the electric 
energy density of the interparticle interactions, i.e. 


a 2 w/ac 2 , may change at some value of the 
concentration. As a result of this feature, which is 
absent in the theory for a dilute case [1], the ability or 
inabilty of the particles to aggregate in the presence of 
applied electric fields depends on the particle 
concentration as well. 

As in the case of conventional atomic systems [81, 
the use of a spinodal curve provides a convenient way 
to distinguish “melastablc” and “unstable” states of a 
suspension subject to an electric field, where we refer 
to the domain between the coexistence curve and the 
spinodal curve in its phase diagram as metastable, and 
the domain beyond the spinodal line as unstable. The 
free energy of a suspension, being considered as a 
function of the particle concentration, is convex in the 
metastable region of the phase diagram “the particle 
concentration-thc electric field strength”, so that the 
free energy will increase with a spontaneous 
concentration fluctuation, giving rise to an energetic 
barrier that will stabilize the suspension. The free 
energy, however, turns into a concave function of the 
particle concentration in the unstable region of the 
phase diagram, so that no energetic barrier to phase 
separation will exist in this state. This distinction can 
correspond to two different mechanisms of clcctric- 
field-induccd transformations in a suspension: spinodal 
decomposition and nuclcation, as occurs in atomic 
systems. For atomic systems, the former (in the 
unstable domain) corresponds to the growth of long- 
wavelength spontaneous concentration fluctuations 
with time whereas the nucleation of microdomains of 
the other phase starts the transformation for the latter 
(in the metastable domain). However, as recent 
research demonstrates (see review in Ref. [91), there is 
no sharp dividing line between nuclcation and spinodal 
transformations in atomic systems. Rather there exists a 
gradual transition in the dynamic behavior of a 
quenched system as the quench point on its phase 
diagram varies from one domain to the other in the 
vicinity of the spinodal curve. 

In any event, we can expect that there exists some 
similarity between how spinodal and nuclcation 
transformations operate in a suspension subject to an 
electric field and how they operate in quenched atomic 
systems, even though the electric-field-induced 
interparticle interactions are anisotropic. In this 
connection, we proceed in Sec. 3 to construct, based on 
Eqs. (7) and (8), a classification scheme which shall 
predict how the topology of the suspension phase 
diagram “the particle concentration-the electric field 
strength” depends on the particle and fluid dielectric 
constants and conductivities and the frequency of 
applied field. This will make it possible to exploit the 
great body of experimental data and theoretical 
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predictions available for the quenching of atomic 
systems and use it as a framework for interpreting the 
morphology and kinetics of aggregation patterns in 
colloidal suspensions generated by the application of 
electric fields. 

3. CLASSIFICATION SCHEME OF DC- 
AND AC-FIELD PHASE DIAGRAMS 

A thorough analysis of Eqs. (7) and (8) which 
yields a complete set of phase diagrams of suspensions 
subject to electric fields as a function of the mismatch 
of the dielectric constants and of conductivities of the 
particles to those of the fluid will be published 
elsewhere. Due to the limitations on the length of the 
paper, we consider only the main features of these 

diagrams shown in Fig. 1 where X = £ f E 2 v p /k B T 

is the electric-to-thermal energy ratio, 1 and 2 are the 
spinodal and coexistence curves, and M and U denote 
the metastable and unstable domains. 



The simplest phase diagram shown in Figs, la 
corresponds to the case when 8 2 W/5c 2 is negative 
over the entire concentration range, 0 < C < C m , so 
that no aggregation at all is predicted to occur in the 
presence of an electric field. In contrast, the phase 
diagram sketched in Figs, lb shows the predicted 

behavior when 3 2 W/3c 2 is positive over 

0 < C < C m so that, regardless of the particle 

concentration, the suspension will start aggregating as 
long as the applied dc field becomes sufficiently large. 
The diagrams depicted in Figs, lc and Id are 

encountered when 3 2 W/0C 2 changes its sign with 
increasing concentration from positive at c = 0 to 
negative at C = C m , so that it equals zero at some 

concentration C r . The phase diagrams given in Figs. 

le and If represent the case when 9 2 W/3c 2 changes 
its sign with decreasing concentration from positive at 
C = C m to negative at C = 0, so that it equals zero at 

some concentration C r . The separation between the 
spinodal and coexistence curves along the vertical lines 
of constant concentrations on the phase diagram in Fig. 

1 increases indefinitely with increasing strength of the 
electric Field. 

As seen from Fig. 1, a metastable domain in the 
high-concentration parts of the phase diagrams in Fig. 

lc (at C r < C < c \ ) and Fig. Id (at C T < C < C m ) as 
well as in the low-concentration parts of the phase 
diagrams in Fig. le (at c\ < C < C r ) and in Fig. If 
(for 0 < C < C r ) is retained regardless of the strength 

of the applied dc field. Hence, only the nuclealion 
mechanisms of electric-field-induced transformations 
will operate in such suspensions under these 
conditions. However, for the phase diagram in Fig. lb, 
an increase in the electric field strength along the 
vertical line of constant concentration will inevitably 
lead to the transition from a metastable to an unstable 
domain at any concentration. Thus for such 
suspensions, the nucleation mechanism will operate in 
weak fields (below the spinodal curve) whereas the 
spinodal decomposition will inevitably operate in 
strong fields (above the spinodal curve). This 
qualitative distinction between the transformation 
mechanisms should manifest itself by the dependence 
of the morphology of the aggregation patterns to be 
formed on the strength of the applied field. 

The map plotted in Fig. 2 demonstrates how the 
topology of the phase diagram of a suspension subject 
to dc fields relates to the mismatch of the dielectric 
constants and the conductivities of the particles and 
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those of the suspending fluid. The domains in Fig. 2 
are lettered to correspond to the phase diagrams in Fig. 
1 whereas the numbers 1, 2, 3, and 4 denote the curves 

a 2 w/a c 2 | c=c =o, a 2 w/a c 2 | c=o = o, 

Jc(a 2 W/0c 2 )dc = O,and 

0 

c m 

J ( c rn “ c)(3 2 w/5c 2 )ic = 0 , respectively, 
o 



Figure 2. 

An increase in the frequency of an ac electric field 
affects the shape and the size of sub-domains in Fig. 2. 
For example, as G)t s increases, the domain “a” located 

in the region (7 p < O f shrinks to the line £ p = E f , so 

that it disappears gradually as CDt s — > 00 . On the 
other hand, the domain “a” located in the region 
O p >O f is retained while its boundaries move 

indefinitely to the right as (0t s — > ©o and, moreover, 
the separation between them along the lines 
£p/e f const increases. In particular, this 


demonstrates why 3 2 W/3c 2 should always become 
positive at any fixed values of the and 

e /e f ratios as the field frequency is made 
sufficiently high. 

It should be pointed out that the particle content, 
C 2 , of the high-concentration phase of a suspension in 
the phase diagrams in Fig. 1 increases dramatically 
with the strength of the applied field. Thus, it is quite 
realistic to suppose that a high-concentration phase in 
the diagrams in Fig. 1 should undergo subsequent 
disorder-to-order transitions in the presence of 
sufficiently strong electric fields as occurs in the 
absence of an electric field [3, 4]. In this connection, 
we also mention recent studies of non-conducting 
suspensions subject to strong electric fields f 1 0] which 
demonstrate that a high-concentrated phase in such 
colloids eventually form a crystalline body-centered 
tetragonal solid. Unfortunately, a description of these 
transformations cannot be included in the diagrams in 
Fig. 1 because of the absence of an equation for the 
free energy of conducting anisotropic aggregates, 
which is required for this analysis. 
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